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ABSTRACT 

With the ultimate aim of distinguishing between various models describing the formation 
of galaxy halos (e.g. radial or multi-phase collapse, random mergers), we have completed a 
spectroscopic study of the globular cluster system of M31. We present the results of deep, 
intermediate-resolution, fibre-optic spectroscopy of several hundred of the M31 globular clusters 
using the Wide Field Fibre Optic Spectrograph (WYFFOS) at the William Herschel Telescope in 
La Palma, Canary Islands. These observations have yielded precise radial velocities (±12 kms -1 ) 
and metallicities (±0.26 dex) for over 200 members of the M31 globular cluster population out 
to a radius of 1.5 degrees from the galaxy center. Many of these clusters have no previous 
published radial velocity or [Fe/H] estimates, and the remainder typically represent significant 
improvements over earlier determinations. We present analyses of the spatial, kinematic and 
metal abundance properties of the M31 globular clusters. We find that the abundance distribution 
of the cluster system is consistent with a bimodal distribution with peaks at [Fe/H] ~ — 1.4 and 
—0.5. The metal-rich clusters demonstrate a centrally concentrated spatial distribution with a 
high rotation amplitude, although this population does not appear significantly flattened and is 
consistent with a bulge population. The metal-poor clusters tend to be less spatially concentrated 
and are also found to have a strong rotation signature. 

Subject headings: galaxy formation, globular cluster systems, M31, spectroscopy, kinematics 

1. Introduction 

The mechanisms involved in galaxy formation 
remain one of the major unsolved problems in as- 
tronomy. Over the past several decades, various 
groups have proposed a range of models which en- 
deavor to explain the observed properties of galax- 
ies and their globular cluster systems (GCSs). One 
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early formation model from Eggen, Lynden-Bcll, 
& Sandage (1962) argues for a single, large-scale 
collapse of material to form galactic bodies such 
as the Milky Way. A principal competing model 
maintains that formation has occurred through 
random mergers of fragmented gas clouds over the 
course of billions of years, implying a hierarchical 
origin (Searle & Zinn 1978; Cole et al. 2000). In or- 
der to interpret the bimodality seen in the globular 
cluster metallicity and color distributions of many 
galaxies (Ashman & Zepf 1998), additional sce- 
narios have been put forth such as multiphase in 
site formation (Forbes, Brodie, & Grillmair 1997), 
major mergers (Zepf & Ashman 1993) and tidal 
stripping/capture (Cote, Marzke, & West 1998). 

Support for the idea that the Milky Way halo 
formed in a hierarchical or episodic manner arises 
from the lack of a significant abundance gradient 
throughout the entire halo globular cluster pop- 



ulation (Parmentier et al. 2000), a notable de- 
crease in the radial component of the velocity el- 
lipsoid of the stellar halo beyond the solar circle 
(Sommer-Larsen et al. 1997), and the detection of 
a possible spread in cluster and stellar ages within 
the outer halo (Sarajedini, Chaboyer, & Demar- 
quc 1997; Rosenberg et al. 1999; Stetson et al. 
1999). In addition, some studies have uncovered 
evidence of kinematic and chemical sub-clustering 
within the globular cluster population of the An- 
dromeda galaxy which may point to accretion rem- 
nants (Ashman & Bird 1993; Saito & lye 2000). 
However, it is possible that some degree of frag- 
mentation might arise in an infall scenario due to 
thermal and gravitational instabilities within the 
collapsing proto-galactic gas cloud. A better un- 
derstanding of the fragment properties including 
chronology, chemistry and, in particular, kinemat- 
ics would help us to identify the origin of sub- 
clustering and to differentiate between the various 
formation models. 

It has been noted often that globular clusters 
(GCs) provide us with ideal probes of galaxy struc- 
ture and formation mechanisms. As the oldest 
galactic stellar systems, globular clusters hold the 
key to uncovering the formation record of their 
host environments. Correlations between chemi- 
cal content (which points to chronology) and the 
kinematics of globular clusters allow us to investi- 
gate the early formative stages of a galaxy and, to 
some extent, probe its ongoing evolution. 

The M31 globular cluster system provides an 
obvious and desirable target of observation: (i) its 
proximity makes it the most accessible large sys- 
tem outside of the Milky Way; (ii) with ;> 435 con- 
firmed candidate members (Barmby et al. 2000), 
the M31 GCS is sufficiently populous that it pro- 
vides a statistically significant sample size; (iii) 
unlike the Milky Way GCS, the M31 population 
can be observed over the full extent of the galaxy 
with less severe effects of line-of-sight contamina- 
tion; (iv) all of the M31 clusters lie at essentially 
the same distance, and thus are not subject to 
such large individual distance uncertainties as are 
the Milky Way globulars. However, the proxim- 
ity of M31 introduces certain issues of contami- 
nation within the cluster sample, such as a diffi- 
culty in distinguishing between GC candidates and 
star-forming regions in many photometric surveys. 
Furthermore, faint globular clusters can easily be 



missed against the bright galaxy background. 

Previous studies of the Milky Way and M31 
globular cluster systems have revealed that these 
two populations exhibit some remarkable simi- 
larities. Zinn (1985) found that the Galactic 
GCS demonstrated clear evidence for bimodality 
in its metallicity distribution, and proposed that 
it could be considered as two distinct subsystems: 
a flattened, metal-rich, rapidly rotating, kinemati- 
cally cold disk system, and a spherical, metal-poor 
halo with lower net rotation and higher velocity 
dispersion. Signs of bimodality have been found 
in the M31 metallicity distribution (Ashman & 
Bird 1993) and in metallicity-sensitive color dis- 
tributions (Barmby et al. 2000). Separating the 
M31 clusters into metallicity sub-systems along 
the same lines as the Milky Way GCS, Huchra, 
Brodie, & Kent (1991) found that the metal-rich 
clusters within R ~ 7' form a rapidly rotating 
(100 — 200 kms -1 ) disk- like system, whereas the 
metal-poor clusters within this radius exhibit no 
net rotation. Beyond this radius, they were unable 
to distinguish the metal-rich from metal-poor pop- 
ulations on the basis of kinematics. Later studies 
have provided additional evidence supporting the 
view that M31 has a rapidly rotating system of 
metal-rich globular clusters (Ashman & Bird 1993; 
Barmby et al. 2000). 

Armandroff (1989) confirmed the presence of a 
modest metallicity gradient within the system of 
disk clusters in the Milky Way. In the first ex- 
tensive spectroscopic survey of the M31 GCS, van 
den Bcrgh (1969) found no significant metallicity 
gradient with position but did report a wide range 
of abundances at any given radius. More recently, 
Huchra et al. (1982, 1991) and Sharov (1988) have 
shown that there is evidence for a weak but mea- 
surable metallicity gradient with projected radius 
in the M31 cluster system. A similar finding was 
reported by Barmby et al. (2000) based on their 
large sample of metal-sensitive colors. Huchra, 
Brodie, & Kent (1991) found a mean M31 GC 
metallicity of [Fe/H] = —1.2, slightly higher than 
that for the Milky Way. 

There remains a lack of persuasive evidence 
that would allow us to distinguish between halo 
formation models — or indeed to suggest a clear 
hybrid of (or alternative to) the classical formation 
scenarios. It is clear that an extensive database of 
high-precision data for the kinematic and chem- 
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ical tracers of a large number of host galaxies is 
required in order to help solve this currently in- 
tractable problem of galaxy formation. Only then 
can the samples be examined for clear evidence 
of kinematic anomalies, mctallicity gradients and 
sub-structure which, if present, may point to the 
formation and enrichment processes which have 
occurred in the galaxy. With this objective, we 
present spectroscopic results for ^ 200 M31 clus- 
ter candidates, adding a large number of new ob- 
servations to the database. 

This paper will take the following form: Sec- 
tion 2 contains a description of the target selection 
and coordinate determinations. In Section 3, we 
provide a summary of the data acquisition and re- 
duction procedure along with a description of the 
velocity and metallicity determinations. We incor- 
porate these new results into the broader sample 
of available M31 cluster data to provide an anal- 
ysis of the abundance distribution of the globular 
cluster system of our nearest large spiral galaxy 
neighbor in Section 4. In Section 5, we examine 
the kinematics of the cluster system and provide a 
mass estimate for the galaxy. A summary and dis- 
cussion of the results of this study are presented 
in Section 6. 

2. Target Selection and Coordinates 

The list of targets was selected from the cat- 
alogue of Battistini et al. (1987), with some ad- 
ditional central targets from unpublished CCD 
photometry by J. Huchra. Cluster positions ac- 
curate to Qf.'2 were determined from Automatic 
Plate Measuring (APM) scans of deep photo- 
graphic plates of a 6° x 6° area around M31. The 
accuracy of the APM cluster coordinates represent 
a marked improvement over existing M31 CCS po- 
sitions, which have typically been measured to the 
nearest arcsecond. 

Galactocentric coordinates for the GCs were 
computed relative to an adopted M31 central po- 
sition of a = 00 /l 40 m 00.1 s , 5 = +40°59'43" 
(B1950): 

X = C\ sin (PA) + C 2 cos (PA) (1) 
Y = -Ci cos (PA) + C 2 sin (PA), 

where C\ = [sin (a — a ) cos 5} and C 2 = [sin 5 cos <5 — 
cos (a — eeo) cos <5 sin Jo] . The X coordinate rep- 
resents the position along the major axis of the 



galaxy increasing towards the north-east along the 
major axis. The Y coordinate is the distance from 
the major axis above and below the galactic plane, 
with positive values towards the north-west. The 
position angle (PA) of the major axis was taken 
to be 38° (Kent 1989). 

Numerous groups have contributed to the iden- 
tification of the globular cluster members of M31 
(Hubble 1932; Vetcsnik 1962; Sargent et al. 1977; 
Crampton et al. 1985; Battistini et al. 1980, 1987, 
1993). The ultimate result of this combined effort 
is that most of the clusters are known by more 
than one name, a fact which has the potential to 
cause confusion. The letters preceding the identi- 
fication number in the M31 cluster target names 
given hereafter indicate the relevant catalogue ref- 
erence as provided in Table 1. 

3. The WYFFOS Data 

The data for this study were obtained at the 
William Hcrschcl 4.2 m telescope 7 (WHT) in 
La Palma, Canary Islands, during the nights of 
November 3 to 6, 1996. The WYFFOS Wide Field 
Fibre Optic Spectrograph (Bridges 1998) was used 
with two gratings to provide a total spectral cov- 
erage of ~ 3700 - 5600 A. The H2400B 2400- 
line grating yielded a dispersion of 0.8 A/pixel 
and a spectral resolution of 2.5 A over the range 
3700 — 4500 A. This wavelength range was se- 
lected to investigate the CN feature at 3883 A, 
the H & K lines of calcium, H6, the CH G-band 
and the 4000 A continuum break. The R1200R 
1200-line grating observations had a dispersion of 
1.5 A/pixel and a resolution of 5.1 A over the 
spectral range 4400 — 5600 A to add absorption 
features such as H/3, the Mgb triplet and two iron 
lines near 5300 A. 

Six different fibre configurations were selected 
to target a total of 288 globular cluster candidates 
over the full extent of the major axis of M31. An 
overlay of the WYFFOS target fields is shown su- 
perimposed on a Digitized Sky Survey image of 
the galaxy in Figure 1. A record of the observa- 
tions is provided in Table 2, which includes the 
coordinates of the field centers, observation dates 

7 The William Hcrschcl Telescope is operated on the island 
of La Palma by the Isaac Newton Group in the Spanish 
Observatorio del Roque de los Muchachos of the Institute 
dc Astrofisica de Canarias. 
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and exposure times. Regrettably, minor axis fields 
were not observed during this run due in part to 
a loss of time from poor weather. 

Spectra of twilight flats, argon arc lamps, and 
velocity template stars were also obtained during 
the same nights of observation as the targets. The 
data reduction — including fibre throughput cor- 
rections, spectral extraction, sky subtraction and 
wavelength calibration — was accomplished using 
the wyf _red WYFFOS Multi-fiber IRAF 8 package 
written by Jim Lewis. Ensemble median RMS val- 
ues for the wavelength calibration were typically 
0.02 - 0.08 A for the H2400B grating images and 
0.08 - 0.13 A for the R1200R grating. 

Most of the globular clusters in the M31 fields 
are superimposed on a relatively bright, non- 
uniform background of galaxy light. For the cen- 
tral WYFFOS fibre configurations, roughly 20 
blank-sky fibres were arranged in a regular pat- 
tern across the 1° field to characterize the back- 
ground intensity. The outer fields contained more 
dedicated sky fibres thanks to a lower density of 
targets in the area. Spectra obtained through 
the dedicated background sky fibres were aver- 
aged to form a composite sky spectrum which was 
then subtracted from each target. Attempting to 
match one or a few nearby sky fibres to each tar- 
get spectrum (i.e. to associate the local sky level 
with each cluster) did not prove advantageous, 
since the resulting losses in signal-to-noise largely 
outweighed any gains in compensating for back- 
ground variability even within the central fields. 
Examples of reduced spectra obtained using both 
gratings are shown in Figures 2 and 3 with the 
absorption features of interest labelled. 

3.1. Radial Velocities 

Radial velocities for the M31 globular clusters 
were calculated by performing standard Fourier 
cross-correlations of each target spectrum with the 
high signal-to-noise spectra of template stars of 
known velocity (Tonry & Davis 1979). Two spec- 
tra were acquired for each target, obtained using 
each of the R1200R and H2400B gratings. Two 
template stars, HD12029 (K2 giant) and HD23169 
(G2 dwarf), were also observed during the run 

Image Reduction and Analysis Facility, distributed by the 
National Optical Astronomical Observatories, which is op- 
erated by AURA under contract with the NSF. 



to yield four cross-correlation combinations for 
the derivation of velocities for each cluster target. 
We required that each target yield a minimum of 
two satisfactory cross-correlations with peak val- 
ues CC > 0.3 to qualify as a successful result (a 
result of CC — 1 would indicate a perfect cor- 
relation). A mean of the cross-correlation results, 
weighted by the associated Tonry & Davis R- value 
and corrected for solar motion, was used to calcu- 
late heliocentric radial velocities for a total of 202 
M31 clusters candidates in the sample. The re- 
mainder represent: 

1. Confirmed non-clusters (foreground stars, 
background galaxies, HII regions, open clus- 
ters etc.) identified by spectral features 
grossly atypical of globular clusters or ra- 
dial velocities that were significantly beyond 
reasonable limits for the M31 GC population 
( | w bs — v sys | > 3<t„ where a v ~ 150 kins" 1 ). 
The status of these objects is noted in the 
final column of Table 3. 

2. Target spectra with insufficient signal-to- 
noise to generate reliable velocities. 

3. Seven objects with inaccurate coordinates 
due to misidentifications or saturation on the 
APM scans. 

The radial velocities obtained from the WYF- 
FOS spectra are provided in Table 3. Over half 
(109) of the M31 GCs for which velocities were ob- 
tained had no previously published spectroscopic 
data. The majority of the remaining clusters 
(72) had published velocities with large uncertain- 
ties of 20 — 80 kins" 1 , and the remainder (21) 
had modest to high-precision velocities with er- 
rors < 20 kms" 1 . The positions of the WYF- 
FOS targets in each of these categories are shown 
schematically in Figure 4. 

Robust fits revealed no significant differences 
between the velocities determined using either 
template object with a given grating. A compar- 
ison of the velocities calculated for a given target 
from the two gratings showed only slightly larger 
differences; the derived velocities found using the 
R1200R grating data, for example, were at worst 
~ 6 kms" 1 larger than those determined from 
the spectra obtained using the H2400B grating. 
An uncertainty of ±12 kms -1 is adopted here for 
the M31 clusters in our sample based on the RMS 
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of the velocity deviations for 28 targets observed 
in multiple fields. On average, this represents an 
overestimate of the true uncertainty since it incor- 
porates certain internal systematic errors as well 
as statistical errors. 

Figure 5 shows a comparison of those ra- 
dial velocities in the present study that overlap 
with other samples from the literature. A lin- 
ear fit, weighted by the inverse square of the 
published velocity uncertainties, yields u pu bi = 
0.43 + 1.02 v bs kms -1 . Overall, there is good 
agreement between the WYFFOS velocities de- 
rived herein and those of other studies within the 
adopted errors. Encouragingly, the RMS residuals 
of independent fits to the separate velocity sam- 
ples were particularly low for the clusters which 
overlap with the high-precision echelle studies of 
Peterson (1989: vp S g = —5.8 + 0.96 w bs kms -1 
with RMS = 5.7 kms -1 for 8 objects) and Dubath 
& Grillmair (1997: v BG = 0.58 + 1.02 v ohs kms" 1 
with RMS = 8.2 kms" 1 for 5 objects). Note that 
the only high-precision result that deviates signif- 
icantly from our fit is for B29, a cluster for which 
Peterson (1989) notes her velocity may be incor- 
rect. This GC was not used in the calculation 
of the weighted fit or RMS, and it is circled in 
Figure 5. 

There are two other points in Figure 5 which 
deviate significantly from the fit, and thus de- 
serve mention. The open diamond above the 
best fit line represents the M31 globular cluster 
B301-S22: Barmby et al. (2000) obtained v r = 
—30 ±20 kms -1 , in poor agreement with the value 
of — 374±12 kms -1 measured in the current study. 
The Barmby et al. coordinates match well with 
the APM coordinates for this object as provided 
in Table 3, so the velocity deviation does not ap- 
pear to be due to a target misidentification in 
either study. This cluster was also measured by 
Federici et al. (1993), who obtained a velocity of 
—419 ± 30 kms -1 , in good agreement with the 
WYFFOS result. The Barmby et al. sample also 
had the largest RMS residuals in the comparison 
fits: RMS = 110 kms -1 for 16 objects, not includ- 
ing B301-S22. The reasons for this large discrep- 
ancy are not known. The open upright triangle 
which deviates from the fit in Figure 5 represents 
B109-S170, for which Huchra, Brodie, & Kent 
(1991) measured a velocity of —613 ± 24 kms -1 
(compared to -372 ± 12 kms -1 from the WYF- 



FOS results). Again, the coordinates for this ob- 
ject match well between both studies and the de- 
viation does not appear to be due to a target 
misidentification. 

3.2. Absorption and Reddening Correc- 
tions 

Each spectrum was corrected for the effects of 
atmospheric absorption and for reddening in both 
the Milky Way and in M31 itself. Color excess 
values upon which the reddening calculations are 
based were determined using the slope (S) of the 
continuum between ~ 4000 — 5000 A as tabulated 
in Crampton et al. (1985), where available: 

E B -v = -0.066 S+ 1.17 (B- V) - 0.32. (2) 

The slope parameter has been shown to be a good 
estimator of reddening, and is not very sensitive 
to cluster metallicity (Elson & Walterbos 1988). 

Targets with published (B — V) colors but with 
no slope parameters or intrinsic color estimates 
were de-reddened to a typical color of (B — V)o = 
0.76 for M31 globular clusters at large galactocen- 
tric radii (Crampton et al. 1985). Clusters with no 
S values (or colors) and those for which the cal- 
culated excess values were lower than a minimum 
foreground color excess of Eb-v = 0.10 (due to 
errors in the slope parameter or in the photom- 
etry) were assigned this minimum value (Frogel, 
Persson, & Cohen 1980; Elson & Walterbos 1988). 

3.3. Abundance Determinations 

To rank the M31 clusters in the WYFFOS sam- 
ple by metallicity, we have measured the strengths 
of various absorption features in their integrated 
spectra. Line indices were calculated to measure 
the signal in a wavelength-delineated region cen- 
tered on each spectral feature relative to the signal 
in the blue and red continuum zones flanking the 
feature. In light of the fact that no flux standard 
star spectra were obtained during the WYFFOS 
observing run, it was not possible to convert pho- 
ton counts directly into flux for these targets. Our 
motivation, however, was to obtain a metallicity 
ranking of the globular clusters in M31; relative 
flux calibration was not essential as its absence 
will not significantly affect narrow features and 
mean count ratios were used instead. Our line 
and continuum bandpasses were defined as shown 
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in Table 4 (Burstein et al. 1984; Brodie & Hanes 
1986; Brodie & Huchra 1990). 

Each target spectrum was shifted to zero ra- 
dial velocity to define consistent line index wave- 
lengths in the rest frame. The velocity disper- 
sion of the M31 globular cluster system is roughly 
150 krns" 1 , a spread which translates into a nom- 
inal range of 4cr v > 600 kms -1 in velocity. This 
would correspond to a wavelength shift of ~ 9A 
at A = 4500 A, a significant fraction of the width 
of most of the bandpasses selected for measure- 
ment. Thus, metallicity determinations were not 
pursued for target spectra which had yielded no 
velocity cross-correlations. 
The generalized feature index is defined as 



-2.5 log 



2Ti 



C2 



(3) 



where Ti , Tc\ and fci represent the mean 
count levels in the feature bandpass and associ- 
ated continuum bandpasses CI and C2 as given by 

T = (A 2 - Ax)- 1 ill TdX ( Brodie & Hanes 1986 ; 
Brodie & Huchra 1990). The continuum break 
feature A = 2.5 log [^02/^01] is a measure of the 
discontinuity in the continuum due to Fraunhofcr 
line-blanketing shortward of 4000 A. This feature 
is somewhat more vulnerable to errors in the red- 
dening and absorption corrections, as it spans a 
broad range in wavelength and is asymmetric. 

Figure 6 shows a comparison of our line in- 
dices with published M31 globular cluster [Fc/H] 
values from various sources (Bonoli et al. 1987; 
Huchra, Brodie, & Kent 1991; Barmby ct al. 2000). 
The RMS differences between the line indices cal- 
culated from spectra observed in more than one 
field through different fibres were used to estimate 
the internal errors for each feature, and these are 
shown in the second column of Table 5. The re- 
sults of linear fits in the form LI = a + 6 [Fc/H] 
weighted by 1 / ^pFc/H] are provided in columns 3 
and 4 of Table 5. Note that only spectroscopic 
metallicities for overlapping targets were consid- 
ered in the calibration and hence the Bonoli ct 
al. (1987) near-infrared photometric metallicities 
were not used in the fits. In column 5 of Table 5, 
we provide the RMS residuals of the fit and in col- 
umn 6 we list the strength of the r coefficient, a 
measure of the linear correlation between the line 
indices and published metallicities. 



To best represent the metal abundances of the 
WYFFOS targets observed, a line index was se- 
lected to contribute to the calibration if it met the 
following four criteria: 

1. The uncertainties associated with the line in- 
dex did not exceed ~ 20% of the full range 
of observed values. The Fe52 index did not 
satisfy this condition and was thus rejected 
from the calibration. 

2. The linear correlation coefficient r exceeded 
a minimum level of r — 0.65. 

3. Numerical experimentation demonstrated 
that the feature strength was stable in 
the presence of significant changes in the 
adopted value of cluster reddening (e.g. in 
comparison with any adopted reddenings 
which differed from unpublished M31 values 
from Barmby et al. (2000)). The continuum 
break was found to be susceptible to redden- 
ing uncertainties and was thus rejected. 

4. Only one line index that characterizes a 
given spectral feature was incorporated in 
the estimate. In the case of the magne- 
sium triplet, the indices MgG, Mg2 and Mgb 
all provided an acceptable measure of es- 
sentially the same absorption feature. Of 
these, the Mgb index demonstrated the high- 
est correlation coefficient and a relatively low 
RMS uncertainty and was therefore selected 
for use over the others. 

Line indices which met the above criteria in- 
clude the CH (G-band) feature, Mgb to charac- 
terize the strength of the magnesium b triplet, and 
the Fe53 iron line. Final cluster metallicities have 
been determined from an unweighted mean of the 
[Fe/H] values calculated from the CH(G), Mgb 
and Fc53 line strengths for 202 of the WYFFOS 
targets. These metallicities and their associated 
errors are provided in columns 7 and 8 of Table 3. 
A plot of the correlation between our [Fe/H] val- 
ues and those in other samples from the literature 
is presented in Figure 7. The best-fit of the data 
in Figure 7 has a slope of 0.94 ±0.02 and the RMS 
of the fit residuals is 0.24. The median value of 
the formal errors on our WYFFOS metallicities is 
±0.26 dex. 
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3.4. The Best Current Sample 

In the subsequent analysis, our results were 
merged with the velocities and [Fe/H] values from 
other sources in order to produce the best available 
sample of spectroscopic data. Wherever data from 
other publications overlapped with our WYFFOS 
sample and demonstrated smaller associated un- 
certainties, these data superseded the WYFFOS 
results; such cases are noted in the final column of 
Table 3. In this final sample, 191 WYFFOS veloci- 
ties and 189 metallicities were combined with data 
from Peterson (1989), Huchra, Brodie, & Kent 
(1991), Federici et al. (1993), Dubath & Grillmair 
(1997), and Barmby et al. (2000) to yield a total 
M31 globular cluster spectroscopic database com- 
prising 321 velocities and 301 metallicities. 

Barmby et al. (2000) identify the following tar- 
gets as having spectral properties which are in- 
consistent with globular clusters: B55-S116 (stel- 
lar), B308 (galaxy), B341-S81 (stellar) and B392- 
S329 (stellar). We concur with the identification 
of B308 as a definite non-cluster and with that of 
B392-S329 as a probable non-cluster, and do not 
carry these through the remainder of the analysis. 
We do not find any clear abnormalities in the spec- 
tra of the other targets and thus cannot at present 
reject these as confirmed non-clusters. 

In an effort to verify that there are no obvi- 
ous biases in the metallicity sample studied here, 
we examine the color distribution of our sample 
and compare it to that of the overall M31 GCS. 
The mean estimated intrinsic color of 412 M31 
globular clusters with (B — V) data from Barmby 
et al. (2000) was (B - V) = 0.67 ± 0.04, while 
that for the sample of those with [Fe/H] values 
used in the metallicity analysis presented herein 
(264 clusters with colors) was determined to be 
(B - V) = 0.65 ± 0.05. Based on this and a 
comparison of the color histograms for the metal- 
licity sample and the overall GCS (Figure 8), we 
conclude that there are no significant population 
biases with respect to the color of the metallicity 
sample presented herein. A Kolmogorov-Smirnov 
two-sample test on the distributions in Figure 8 
confirms what the eye suggests: the color distri- 
bution of the metallicity sample is drawn from the 
same parent population as the overall color sample 
at a high confidence level (99.5%). 

Clearly, the globular clusters at small galacto- 



centric radii will suffer from more severe effects of 
incompleteness due to obscuration by the disk and 
bulge of M31. Depending on the color/metallicity 
distribution as a function of position within the 
galaxy, this observational effect will conceivably 
bias the observed metallicity sample. We will re- 
turn to this point in Section 4.1. 

4. M31 Cluster Metallicities 

The metal abundance properties of M31's glob- 
ular clusters are of interest as they bespeak the 
galaxy's formative processes and enrichment his- 
tory. The metallicity distribution for the M31 
GCS is presented in Figure 9, along with a simi- 
lar plot for the Milky Way GCS for comparison. 
The mean metallicity of the full complement of 301 
M31 globular clusters is [Fe/H] = -1.21 ± 0.03, 
comparable to the value of —1.27 ± 0.05 obtained 
for the Milky Way GCS (based on data sup- 
plied in Harris (1996), revised 1999). However, a 
Kolmogorov-Smirnov test reveals that the Galac- 
tic and M31 globular clusters have [Fe/H] distri- 
butions which are drawn from the same parent 
population at only the 55% confidence level. 

As is frequently the case with binned data, the 
appearance of the histogram can be ambiguous 
and potentially misleading, and thus a more ro- 
bust method of analysis is desirable. In order to 
examine the shape of the metallicity distribution 
without relying on binning methods, we turned 
to the KMM algorithm used by Ashman, Bird, & 
Zcpf (1994). KMM mixture modeling operates un- 
der the assumption that the data is independently 
drawn from a parent population which consists of 
a mixture of N Gaussian distributions. The KMM 
test has previously been applied to the M31 glob- 
ular cluster metallicity distribution by Ashman 
& Bird (1993), albeit with a significantly smaller 
data sample. 

A bimodal distribution was tested against the 
assumption of a unimodal shape (the null hypoth- 
esis) for the GCS metallicities. Since there is no 
obvious reason to assume that both [Fe/H] groups 
in a bimodal fit would exhibit a common co- 
variance, hctcroscedastic fitting was applied such 
that the covariances for the two distributions were 
not constrained to be identical. The KMM re- 
sults revealed that this 2-group fit demonstrated 
a marginal improvement over a single group fit, 
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although Ashman, Bird, & Zepf (1994) cautioned 
that the output probability in the heteroscedastic 
case is more difficult to interpret and may, in and 
of itself, be suspect. We repeated the KMM test 
on more than 1500 simulations with bootstrap re- 
sampling and found that the bimodal case was pre- 
ferred over the unimodal distribution at a median 
confidence level of higher than 97%. The bimodal 
test returned Gaussian fits characterized by means 
at [Fe/H] = -1.44 (a 2 = 0.22) for the metal-poor 
distribution and at [Fe/H] = -0.50 (a 2 = 0.13) 
for the metal-rich clusters. The a posteriori prob- 
abilities of group membership returned by KMM 
assigned 231 clusters to the metal-poor distribu- 
tion and 70 to the metal-rich population. 

In their analysis, Elson & Walterbos (1988) 
treated the M31 GCS as comprising three distinct 
metallicity sub-populations: metal-poor, interme- 
diate, and metal-rich. We employed the KMM test 
to ascertain whether a trimodal distribution was 
also consistent with the measured [Fe/H] values 
for the M31 cluster system. The trimodal KMM 
results with bootstrap resampling demonstrated 
that a 3-group fit was indeed preferred over a sin- 
gle Gaussian distribution at a relatively high con- 
fidence level (96%). The output population means 
obtained for the best 3-group fit were located at 
[Fe/H] = -1.59 (a 2 = 0.19), -1.18 (a 2 = 0.18) 
and -0.47 (cr 2 = 0.12). The numbers of clusters 
assigned to each group were N =153, 76 and 72, 
respectively. 

In light of the variation in the magnitude of 
[Fe/H] uncertainties quoted within the WYFFOS 
sample and others, it is appropriate to test the ef- 
fects of these uncertainties on the inferred shape of 
the distribution. There are 229 GCs with metal- 
licity uncertainties lower than 0.45 dex, a value 
roughly equal to the variance of one of the pop- 
ulations as determined from the bimodal results. 
We applied the KMM test to this higher-precision 
sample to test the stability of our earlier conclu- 
sions. The bimodal case was again found to rep- 
resent a statistically significant improvement over 
the unimodal assumption with peaks similar to 
those found when we incorporated the complete 
sample. The input assumption of three peaks, 
however, did not return adequate fits: the a pos- 
teriori probabilities of group membership in the 
trimodal case allocated only two clusters to the 
intermediate population. 



There is an unavoidable smearing of the metal- 
licity distribution due to metallicity calibra- 
tion uncertainties, absorption and reddening of 
bulge/disk clusters and other line-of-sight prob- 
lems, or errors incurred as a result of combining 
different datasets. As a further test, simulations 
were run in which we shuffled the high-precision 
[Fe/H] sample by a random fraction of their as- 
sociated uncertainties. This returned a result of 
marginal significance (median p- value ~ 0.1) for 
the 2-group case, and an insignificant result (p- 
value ~ 0.4) for the 3-group heteroscedastic trial. 
Note that the shapes of the 2-group Gaussian fits 
were once again consistent with those of the orig- 
inal bimodal test using the full sample. 

In light of these KMM results, we tentatively 
reject the possibility of a trimodal GC metallicity 
distribution for the purposes of the analysis which 
follows. While it may yet be reasonable to sepa- 
rate the cluster system into more than two pop- 
ulations — say, by culling out the intermediate 
targets with group membership probabilities be- 
tween 0.4 and 0.6 in the bimodal case — we elect 
not to do so here, as it is not clear that these rep- 
resent a bona fide intermediate population. Based 
on the available data, we argue that the metallic- 
ity distribution of the M31 globular cluster system 
is not unimodal and better resembles a bimodal 
distribution similar in general shape to that of 
the Milky Way (Zinn 1985). Therefore, "metal- 
poor" and "metal-rich" populations shall refer to 
those objects separated into two groups based on 
the a posteriori probabilities returned from the 2- 
group heteroscedastic KMM test with the com- 
plete available sample. Each population was inde- 
pendently fit using an adaptive kernel smoothing 
procedure 9 and the results are shown in Figure 10. 

The positions of our adopted metallicity peaks 
([Fe/H] = -1.44 and -0.50) differ only slightly 
from those quoted by Barmby et al. (2000), who 
found a bimodal distribution in M31 GC metallic- 
ities and mctallicity-sensitive colors with peaks at 
[Fe/H] = —1.4 and —0.6. Our peaks are also sim- 
ilar to the values of [Fe/H] = —1.5 and —0.6 ob- 
tained by Ashman & Bird (1993) using the KMM 
test on an earlier sample of 144 M31 globular clus- 
ters; their study yielded a ratio of metal-rich to 

9 Adaptive kernel smoothing code kindly provided by Karl 
Gebhardt. 
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metal-poor globular clusters of 49:95, in compar- 
ison with a ratio of 70:231 found in the current 
study. 

4.1. Spatial Distribution 

The projected spatial distributions of the 
metal-poor and metal-rich cluster populations are 
shown in Figure 11. The metal-rich population 
is clearly more centrally concentrated, an effect 
which has been demonstrated previously with a 
smaller dataset (Huchra, Brodie, & Kent 1991). 
In the Milky Way, the metal-rich clusters have 
been ascribed as belonging to a thick disk system 
(Zinn 1985; Elson & Walterbos 1988; Armandroff 
1989). Later work has suggested that this is better 
described as a bulge/bar cluster system (Minniti 
1995; Cote 1999; Forbes, Brodie & Larsen 2001). 
The spatial concentration and distribution of the 
metal-rich clusters in M31 appear to be consistent 
with the latter interpretation. 

There is a notable deficit of metal-poor clusters 
in the innermost radial bins shown in Figure 12. 
The already strong central concentration in the 
metal-rich population might mask any analogous 
deficit in this population. A similar lack of cen- 
tral clusters was reported by Elson & Walterbos 
(1988) within the inner 15' of the GCS and is at 
least partially attributable to sampling restrictions 
towards the middle of the galaxy, as well as to ob- 
servational biases incurred as a result of obscura- 
tion by the bulge and disk. The cluster sample at 
large galactocentric radii is bound to be more com- 
plete at faint magnitudes due to the lower degree 
of absorption as compared with the inner cluster 
system. 

Examination of Figure 13 reveals a modest lack 
of clusters fainter than V w 17 within the inner 
10' or so of the GCS. We attempt to estimate the 
incompleteness at small radii in the following man- 
ner. Using available V photometry (Barmby et al. 
2000), we plot the observed luminosity functions 
(LFs) for the metal-rich and metal-poor clusters 
beyond a galactocentric radius of 15'. These outer 
LFs are normalized to the same scale as the dis- 
tributions within the inner 15' by summing under 
the brighter end of the inner and outer histograms 
(V < 16.75 mag) and multiplying by the ratio of 
these values. The observed histograms for the in- 
ner metal-rich and metal-poor clusters are shown 
by the solid lines in Figure 14, and the scaled 



outer histograms are represented by the dotted 
lines. Note that this method provides only a very 
crude estimate of the number of missing clusters at 
V > 16.75 as it presumes, perhaps falsely, that the 
shape of the globular cluster LF does not change 
within the inner region of the galaxy. Dynami- 
cal destruction will have a more severe effect on 
the inner clusters and will thus change the over- 
all shape of the luminosity function. Evidence for 
this effect has previously been noted by Kavelaars 
& Hanes (1997). 

From the scaled outer LFs in Figure 14, we in- 
fer that there are possibly <~ 30 missing metal- 
poor clusters at V > 16.75, whereas there may be 
only ~ 8 faint, metal-rich clusters which are miss- 
ing within the inner 15' of the galaxy. Therefore, 
sample incompleteness is a reasonable explanation 
for the dip in the central bins of the metal-poor 
radial profile. 

Comparing simple counts of the cluster popula- 
tions within quadrants centered on the major and 
minor axes reveals that there is a slight enhance- 
ment in the minor axis fields for the metal-rich 
population: 40 clusters in the East- West sectors 
versus 30 in the North-South. If one takes into ac- 
count only the inner 30' of the GCS, this enhance- 
ment largely vanishes for the metal-rich clusters: 
28 clusters are found within the major axis sectors 
compared to 26 along the minor axis. Based on 
their photometric metallicity estimates, Elson & 
Walterbos (1988) discovered that their sample of 
clusters classified as metal-rich defines a flattened 
system with an ellipticity of e J> 0.4 within a major 
axis radius of 30'. Our cluster counts are not con- 
sistent with such a large degree of flattening of the 
metal-rich system in this region. At present, we 
cannot provide reasonable limits on the flattening 
of the metal-poor cluster distribution, since the 
arrangement of the WYFFOS fields has resulted 
in a selection bias in favor of clusters in major axis 
sectors at radii beyond ~ 30'. 

4.2. Metallicity Gradient 

Several previous works have suggested the pres- 
ence of a weak metallicity gradient in the M31 
cluster system out to radii of ~ 10 kpc (Sharov 
1988; Huchra, Brodie, & Kent 1991). In Fig- 
ure 15, we plot the metallicity of the M31 globular 
clusters as a function of galactocentric radius for 
the full sample and the separate metallicity sub- 
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populations. The large scatter in these plots is due 
at least in part to the fact that we can only mea- 
sure projected positions, but nonetheless we find 
that the mean [Fc/H] does indeed decrease with 
radius. 

A sliding bin method was used to fit the mean 
metallicity in 12' bins centered on each data point, 
with results shown in Figure 15. Robust estimates 
of the mean metallicity in radial bins are consis- 
tent with the presence of a modest gradient in 
[Fe/H] within the full sample. Separating the data 
into sub-populations reveals no clear evidence for 
a gradient in the metal-rich GCS, although the 
spread in metallicity does seem to increase to- 
wards the center of the galaxy. There is an ap- 
parent trend towards decreasing metallicity within 
both the full sample and the inner metal-poor clus- 
ter system, with gradients of roughly —0.017 and 
—0.015 dex/arcmin, respectively. Relatively poor 
sampling of the clusters at large radii makes it 
difficult to comment decisively on the results be- 
yond a radius of about 60', although the metal- 
poor sample appears to level off at a lower mean 
metallicity beyond this distance. 

5. GCS Kinematics 

The maximum likelihood value for the mean 
cluster radial velocity is (v) = —284 ± 9 kms -1 , 
somewhat higher than M31's systemic velocity 
of -300 ± 4 kms" 1 (RC3). The globular clus- 
ters have an overall velocity dispersion of a v = 
156 ± 6 kms -1 . The kinematics of metal-rich 
and metal-poor GCS populations will be discussed 
later in this section. First, we examine the global 
kinematic properties of the M31 GCS. 

5.1. GCS Rotation 

A histogram of the cluster line-of-sight veloci- 
ties is presented in Figure 16. There is strong ev- 
idence for bimodality in the velocity distribution, 
a feature which is attributable to the fact that the 
cluster system is rotating with respect to our line 
of sight. Considering, for the moment, a purely ro- 
tating disk population of globular clusters in the 
reference frame of the galaxy: 

vgg = "rot sin (i) sin (6 + 6) + v„ (4) 
= v' + v G 



where w ro t is the rotation velocity of the popula- 
tion, i is the inclination of the rotation axis to the 
line of sight (77.7°), 9 is the angle about the angu- 
lar momentum axis, 8 is an arbitrary phase offset, 
and v a incorporates the velocity dispersion plus 
random measurement errors. The combined prob- 
ability distribution P(«gc) f° r the M31 globular 
cluster radial velocities is the convolution of P(v') 
with P(v a ). The probability distribution P(v a ) is 
Gaussian and P(0) is uniform, thus P(v') is given 
by: 

PW) = P{e)— (5) 
i i 

27r /~2 . 2 ■ ~~a 

V<t sm * - v 

This represents a bimodal distribution with max- 
ima at ±w rot sin (i) about the systemic velocity 
of the population. The peaks of the (unshiftcd) 
velocity histogram in Figure 16 lie at —415 and 
— 170 kins" 1 , yielding an estimated rotation am- 
plitude of v ro t ~ 125 kms -1 . 

The fact that M31 is rotating has been estab- 
lished by various studies (see Hodge 1992, and ref- 
erences therein). Using a sample of ~ 150 M31 
clusters as dynamical tracers, Huchra, Brodie, & 
Kent (1991) noted a rapid rotation of the hand- 
ful of inner (R <J 7') clusters in their sample 
with [Fe/H] > —0.8 (see also Huchra, Stauffer, & 
Van Speybroeck (1982)). Huchra, Brodie, & Kent 
(1991) discovered that the metal-poor population 
had no rotation signature within the inner galaxy, 
but shared a moderate rotation of ~ 60 kms -1 
with the bulk of the GCS beyond this radius. With 
our increased sample of 321 velocities, we can now 
provide a more comprehensive look at the kine- 
matics of the M31 globular cluster system. 

Perhaps not surprisingly, there appears to be 
no significant rotation of the GCS about the ma- 
jor axis of M31 as shown in the left panel of Fig- 
ure 17. A robust fit to the targets within the in- 
ner 50' along the minor axis (dropping 11 objects 
at large radii) gives v r = —286 + 0.13 Y kms -1 , 
where Y is the distance along the minor axis in ar- 
cminutes. Such a small deviation from a slope of 
zero may simply be due to statistical imprecision 
or variation in the position angle for the galaxy 
(Hodge & Kcnnicutt 1982). 

A plot of the observed cluster radial velocities 
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as a function of projected distance along the ma- 
jor axis of M31 is shown in the right panel of Fig- 
ure 17. A basic 2a culled robust linear fit to the 
velocities along the major axis reveals that the 
globular cluster system indeed has a significant 
rotational component about the minor axis. This 
rotation will now be investigated in more detail. 

We used an adaptive binning technique to gen- 
erate the GCS rotation curve shown in Figure 18. 
Before fitting the rotation curve, the globular clus- 
ter coordinates were transformed into the refer- 
ence frame of the galaxy and the data were folded 
about each axis. The function v = v sys + v Iot sine/) 
was fit to the radial velocities of the clusters 
within radial bins outwards from the galaxy cen- 
ter (Kissler-Patig & Gcbhardt 1998). The angle 
</> represents the position angle of the target ob- 
ject in the projected frame of the galaxy, and v sys 
is the mean velocity of the cluster system. The 
kinematic properties were determined in a series 
of annuli centered on each cluster and encompass- 
ing a specified number of data points (79 points, 
with a minimum of 39 at the extrema). With 
this method, we obtain a rotation amplitude of 
136 kms -1 where the curve flattens off beyond 
R = 30'. 

We then employed the robust biweight estima- 
tors of location and scale (Beers, Flynn, & Geb- 
hardt 1990; Hoaglin, Mosteller, & Tukey 2000) 
as measures of mean velocity and dispersion in 
10 arcminute bins along the major axis of the 
galaxy. The biweight results and their 90% con- 
fidence intervals (based on bootstrap errors) are 
presented in Figure 19 for the full sample and the 
two metallicity sub-populations. Based on the ap- 
pearance of these plots, it seems reasonable to ap- 
proximate the shape of the curve as consisting of 
solid-body rotation within the central 25', with 
a flattening beyond this radius. One-half of the 
velocity difference between the flat regions of the 
curve in Figure 19 yields a rotation amplitude of 
v I0 t — 138 ± 13 kms -1 for the full sample. 

Significant differences between the kinematic 
properties of the metal-rich and metal-poor clus- 
ters may indicate that these populations under- 
went dissimilar formation mechanisms. The max- 
imum likelihood value for the mean velocity of 
the 70 clusters in the metal-rich sample is (v) = 
— 260 ± 18 kms -1 with a dispersion of a v = 146 ± 
12 kms -1 . For comparison, the larger metal-poor 



sample of 231 clusters demonstrates a somewhat 
lower mean velocity of (v) — —290 ± 10 kms -1 
and a slightly higher dispersion of u v = 155 ± 
7 kms -1 , although these values are still consis- 
tent with those of the higher metallicity sample 
within the quoted errors. From the fits of the bi- 
weight results, we determine a rotation amplitude 
of v Iot = 160 ± 19 kms -1 for the metal-rich pop- 
ulation. This value is somewhat larger than the 
metal-poor amplitude of v ro t = 131 ± 13 kms -1 , 
although these values are also not inconsistent 
within their formal errors. 

For each cluster, the projected position angle 
was determined as the angle of maximum positive 
rotation obtained from the maximum likelihood si- 
nusoid fits (Kissler-Patig & Gebhardt 1998). Plots 
of the variation of position angle with radius for 
the full sample and metallicity sub-populations arc 
presented in Figure 20. The input GC positions 
already incorporate a de-rotation of PAjvi3i = 38° 
for the galaxy and, all else being equal, one would 
expect to find that the PA for the cluster system 
is roughly zero. Despite the sizeable spread at 
low radius in the full sample, there is a relatively 
clear increase in PA within the central region of 
the metal-poor cluster system. This effect is con- 
sistent with studies of the central bulge and inner 
disk region of M31 which demonstrate a similar 
increase in PA (Hodge & Kennicutt 1982; Kent 
1989). The M31 GCS shows a notable decrease 
in the PA of the rotation axis beyond about 30'. 
Interestingly, the metal-rich clusters seem to be 
rotating about an axis that is tilted at a slightly 
lower PA with respect to the major axis through- 
out the full radial extent of the population. 

5.2. The Mass of M31 

The kinematics of M31's GCS provide us with 
a dynamical probe of the underlying mass distri- 
bution of the host galaxy (Federici et al. 1990, 
1993). The projected mass estimator (PME) was 
used to place constraints on the underlying mass 
of the system (Bahcall & Tremaine 1981; Heisler, 
Tremaine, & Bahcall 1985): 

N 

M PM e = ^$>^, ( fi ) 

i=\ 

where N is the number of bodies, G is the gravi- 
tational constant, r is the projected galactocentric 
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radius and v is the radial velocity of the body rel- 
ative to the systemic velocity of the galaxy. The 
factor f p depends on the assumptions made re- 
garding the distribution of orbits, and is equal to 
32/-7T and 64/-7T for isotropic and radial orbits, re- 
spectively. 

Admittedly, this projected mass estimate is 
hampered by our limited sample size and coverage, 
as well as a lack of detailed understanding of the 
cluster orbits. Recent work by Evans & Wilkinson 
(2000) has shown that the dynamical tracers in the 
M31 halo exhibit a predominantly isotropic veloc- 
ity dispersion. Applying the PME to the available 
cluster data and assuming isotropic orbits within 
an extended mass distribution, we estimate a mass 
of M M3 i = (4.1 ± 0.1) x 10 11 M using 319 dy- 
namical targets out to a radius of ~ 120' from the 
galaxy center (~ 27 kpc at a distance of 770 kpc). 
This result is comparable to previous mass esti- 
mates obtained out to similar radii as summa- 
rized in Table 6 of Evans & Wilkinson (2000), and 
to the mass of M M3 i ~ (3.1 ± 0.5) x 10 11 M 
found by Huchra, Brodie, & Kent (1991) based on 
their sample of 150 globular clusters. In the (un- 
likely) case of purely radial orbits, our mass esti- 
mate would be increased to twice the value quoted 
above. 

6. Discussion and Conclusions 

We have presented a sample of new spec- 
troscopic observations of M31 globular clusters, 
contributing <~ 200 high-precision velocities and 
[Fe/H] values for what represents a large fraction 
of the dynamical tracers of the galaxy. A compre- 
hensive analysis of the kinematics and abundance 
properties of M31's GCS demonstrates some inter- 
esting similarities and differences when compared 
with the Milky Way's globular cluster population. 

The Galactic and M31 globular cluster systems 
have comparable mean metallicities of [Fe/H] ~ 
— 1.2 to —1.3, although a Kolmogorov-Smirnov 
test provides no conclusive evidence that the two 
distributions are drawn from the same parent pop- 
ulation. A KMM mixture-modeling test reveals 
that the M31 GCS metallicity distribution can 
convincingly be separated into two groups, a bi- 
modality which is analogous to that found in the 
Milky Way GCS. The peaks in the M31 metallic- 
ity distribution lie at [Fe/H] = -1.44 and -0.50. 



A KMM test on the Milky Way GCS metallici- 
ties from the McMaster catalogue (Harris 1996) 
returns peaks at [Fe/H] = —1.59 and —0.56, in- 
dicating that the metal-poor Galactic globulars 
have, on average, a somewhat lower metal abun- 
dance than their M31 counterparts. Of the 301 
M31 globular clusters with available metallicities, 
KMM assigns 70 to the metal-rich peak. This 
represents a somewhat smaller fraction of high- 
metallicity clusters than is found in the Milky Way 
system, for which 45 of 145 are identified with the 
metal-rich population. 

M31's metal-rich globular clusters appear to 
constitute a distinct kinematic subsystem that is 
quite spatially concentrated, consistent with mem- 
bership in a bulge population akin to the in- 
ner metal-rich clusters of the Milky Way (Minniti 
1995). Furthermore, the velocity dispersion of the 
metal-rich cluster system, a v = 146 ±12 kms^ 1 , is 
indistinguishable from the mean bulge dispersion 
of <j v (bulge) = 146 ±6 kms -1 as determined from 
the kinematics of faint planetary nebulae (Lawrie 
1983; van den Bergh 1991). The mean metal- 
licity of this population is quite similar to their 
Milky Way counterparts and demonstrates no ob- 
vious abundance gradient with galactocentric ra- 
dius. Contrary to observations in the Milky Way, 
however, we find no clear signs of flattening in the 
metal-rich component of the M31 GCS despite our 
relative oversampling of major axis fields in the 
current spectroscopic sample. We find that the 
rotation axis of the metal-rich cluster population 
is tilted at a slightly lower position angle: it is 
offset by about 5 — 10° with respect to that of 
the minor axis of M31 and the bulk of its GCS. 
This observation may point to a disparate origin 
for these clusters. Further searches for potential 
cluster candidates hidden by the bulge of M31 will 
be helpful in providing a clearer understanding of 
the kinematics and nature of this inner system. 

In M31, there are more than three times the 
number of metal-poor clusters as there are metal- 
rich. It is not clear if it is possible (or even nec- 
essary) to further separate this larger, more spa- 
tially distended metal-poor population into sub- 
components. If the metal-poor cluster distribution 
does indeed incorporate a thick disk component 
in addition to halo clusters, this would account 
for the relatively high net rotation amplitude 
and large velocity dispersion observed therein. 
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While grouping the clusters based on mctallic- 
ity remains an instructive exercise, [Fe/H] alone 
is not an adequate discriminator of membership 
in the spatially- or kinematically-dcfined compo- 
nents of the galaxy. This is primarily a result of 
the substantial spread and measurable gradient in 
[Fe/H] within the metal-poor cluster system. The 
task of assigning globular clusters to their appro- 
priate component (e.g. bulge, thin or thick disk, 
inner or outer halo) remains somewhat simpler 
within the Milky Way GCS, since here we have 
three-dimensional spatial information rather than 
projected positions. 

We find evidence for a radial metallicity gradi- 
ent (—0.015 dex/arcmin) in the metal-poor popu- 
lation of M31 out to about 60' from the galaxy cen- 
ter. Despite a notable lack of cluster data beyond 
this radius, it does seem that mean metallicity lev- 
els off, although the scatter remains large. This 
gradient in the metal-poor clusters is not incon- 
sistent with a single collapse scenario, yet we rec- 
ognize that a metallicity gradient, in and of itself, 
is not sufficient to confirm an early dissipational 
collapse. Such a gradient could be mimicked by a 
hierarchical halo formation if the metallicity of the 
accreting fragments correlates with their mass and 
mean density (Freeman 1996). The search for sub- 
structure in position/velocity/metallicity parame- 
ter space is clearly of interest to investigate the 
accretion history of the galaxy. An analysis of 
potential sub-clustering within the M31 GCS will 
follow in a subsequent paper. 

A comprehensive investigation of the kinemat- 
ics of the dynamical tracers in M31 would not be 
complete without the addition of planetary nebu- 
lae data. Fortunately, a study of a large sample 
of M31 is currently underway (Halliday et al., in 
preparation. Moreover, these data can be com- 
bined with the globular cluster sample to bet- 
ter probe the underlying mass distribution of the 
galaxy (Wilkinson et al., in preparation). An 
analysis of M31 cluster ages is also in progress 
(T. Bridges, private communication) which should 
provide a better sense of the chronology of GC 
formation within the galaxy. Despite its long and 
rich history of investigation, there is little doubt 
that there are missing globulars in M31's database 
which are awaiting discovery and identification. 
Continuing searches and surveys for M31 cluster 
candidates (Lee et al. 2001) promise to help rem- 



edy this situation and to further fuel this "boom- 
ing industry". It is clear that detailed studies of 
the kinematics, metallicities and ages of the glob- 
ular cluster systems of a large sample of galaxies 
can provide us with the clues necessary for a more 
complete understanding of the mechanisms which 
form and evolve galaxies of all kinds. 
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Table 1 

Target catalogue references 



Code 


Reference 


B 


Battistini ct al. (1980, 1987) 


S 


Sargent et al. (1977) 


BA 


Baade & Arp (1964) 


BoD 


Battistini et al. (1987) Table VI 


DAO 


Crampton ct al. (1985) 


NB 


Battistini et al. (1993) 


V 


Vetesnik (1962) 



Table 2 
Observing log 



Field ID 


aB1950 


&1950 


Night 


UT 


Exp (Sec) 


Grating 


Central 1 


00:40:00.18 


+40:59:59.9 


Nov. 


3/4 


21:00 


6 


X 


1800 


H2400B 


Central 1 


00:40:00.43 


+41:00:00.5 


Nov. 


3/4 


00:43 


3 


X 


1200 


R1200R 


Northeast 1 


00:42:29.88 


+41:36:02.0 


Nov. 


4/5 


21:52 


6 


X 


1800 


H2400B 


Northeast 1 


00:42:29.88 


+41:36:01.1 


Nov. 


4/5 


01:22 


3 


X 


1200 


R1200R 


Southwest 1 


00:37:30.15 


+40:24:01.7 


Nov. 


5/6 


19:32 


3 


X 


1200 


R1200R 


Southwest 1 


00:37:30.51 


+40:24:00.0 


Nov. 


5/6 


20:42 


4 


X 


1800 


H2400B 


Northeast 2 


00:44:59.97 


+42:11:53.4 


Nov. 


5/6 


23:43 


5 


X 


1800 


H2400B 


Northeast 2 


00:45:00.70 


+42:11:56.3 


Nov. 


5/6 


02:40 


3 


X 


1200 


R1200R 


Central 2 


00:39:59.20 


+41:00:05.2 


Nov. 


6/7 


19:26 


3 


X 


1200 


R1200R 


Central 2 


00:39:59.34 


+41:00:03.7 


Nov. 


6/7 


20:38 


4 


X 


1800 


H2400B 


Southwest 2 


00:34:59.80 


+39:47:54.6 


Nov. 


6/7 


23:33 


5 


X 


1800 


H2400B 


Southwest 2 


00:35:00.56 


+39:47:54.6 


Nov. 


6/7 


02:25 


4 x 


1200 


R1200R 



Table 3 
WYFFOS RESULTS 



GC RA (B1950) Dec (B1950) X (') Y (') v (kms" 1 ) [Fc/H] Notes 
(Complete table available in published version) 
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Table 4 
Line Indices and Colors 



Index CI I C2 



HK 


3910 


.00- 


3925 


,00 


3925 


,00 - 


3995, 


.00 


3995 


.00- 


4015.00 


CNR 


4082. 


.00 - 


4118. 


.50 


4144. 


,00 - 


4177. 


.50 


4246 


.00- 


4284.75 


Cal 


4200 


.00 - 


4215, 


.00 


4215, 


,00 - 


4245. 


.00 


4245 


.00- 


4260.00 


CH(G) 


4268 


.25 - 


4283. 


.25 


4283, 


.25 - 


4317. 


00 


4320 


.75 - 


4335.75 


A 


3800 


.00 - 


4000. 


.00 










4000. 


.00- 


4200.00 


H0 


4829. 


.50- 


4848. 


.25 


4849, 


.50 - 


4877. 


.00 


4878 


.25 - 


4892.00 


MgH 


4897 


.00- 


4958. 


,25 


5071, 


.00 - 


5134. 


.75 


5303 


.00- 


5366.75 


Mg2 


4897 


.00- 


4958. 


,25 


5156. 


.00 - 


5197. 


.25 


5303 


.00- 


5366.75 


MgG 


5125 


.00- 


5150. 


,00 


5150, 


.00 - 


5195. 


.00 


5195 


.00- 


5220.00 


Mgb 


5144 


.50- 


5162. 


,00 


5162, 


.00 - 


5193. 


.25 


5193 


.25 - 


5207.00 


Fe52 


5235 


.50 - 


5249. 


,25 


5248, 


.00 - 


5286. 


.75 


5288. 


.00- 


5319.25 


Fe53 


5307 


.25 - 


5317. 


.25 


5314, 


.75 - 


5353. 


50 


5356 


.00- 


5364.75 



Table 5 

Correlation between line indices and [Fe/H] 



Line 


a 


a 


b 


RMS 


r 


H&K 


0.10 


0.51 


0.16 


0.12 


0.37 


CNR 


0.06 


0.18 


0.11 


0.08 


0.51 


Cal 


0.06 


0.00 


-0.01 


0.06 


-0.05 


CH(G)* 


0.08 


0.25 


0.09 


0.06 


0.68 


A 


0.20 


1.31 


0.24 


0.14 


0.70 


H/3 


0.02 


0.03 


-0.05 


0.05 


-0.34 


MgH 


0.04 


0.13 


0.02 


0.04 


0.52 


Mg2 


0.05 


0.20 


0.11 


0.05 


0.72 


MgG 


0.03 


0.13 


0.06 


0.02 


0.75 


Mgb* 


0.04 


0.20 


0.09 


0.03 


0.80 


Fe52 


0.04 


0.11 


0.05 


0.02 


0.61 


Fe53* 


0.03 


0.09 


0.04 


0.02 


0.66 



Note. — Features marked with an asterisk 
were used in the final [Fe/H] calibration. 
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h 50 m 45 m 40 m 35 m 30 m 

Right Ascension (B1950) 

Fig. 1. — Five 1° WYFFOS fields shown superimposed over an optical image of the galaxy from the Digitized 
Sky Survey. Two central fibre configurations were used to yield spectral observations for a total of 288 cluster 
candidates. North is to the top, east is towards the left. 
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3800 3900 4000 4100 4200 4300 4400 



Wavelength (Angstroms) 

Fig. 2. — A sample reduced spectrum for M31 globular cluster B19, obtained using the WYFFOS system 
and the H2400B grating. The ordinate axis is linear in arbitrary units (zero counts is at the bottom). 




4250 
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5000 



5250 



5500 
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Fig. 3. — A sample reduced spectrum for M31 globular cluster B19 obtained using the WYFFOS system 
and the R1200R grating. The ordinate axis is linear in arbitrary units (zero counts is at the bottom). 
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Fig. 4. — The positions of the 202 M31 GCs for which radial velocities were obtained from the WYFFOS 
spectra. Over half of this sample had no previously published spectroscopic data. Of those that had, our 
velocities (with adopted uncertainties of ±12 kms -1 ) represent significant improvements over most of the 
pre-existing spectroscopic data. 
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WYFFOS Radial Velocity (kin s l ) 



Fig. 5. — A comparison between radial velocity results obtained from this study and those in other publi- 
cations: P89 = Peterson (1989); DG = Dubath & Grillmair (1997); BOO = Barmby ct al. (2000); HBK = 
Huchra, Brodie, & Kent (1991); HSV = Huchra, Stauffer, & Van Speybroeck (1982); F93 = Federici et al. 
(1993). Filled points represent high-precision velocity determinations with quoted errors less than 20 kms -1 . 
The unity relation is represented by the dotted line, and the results of a linear fit with a 1 /a 2 weighting is 
shown by the solid line. The circled (filled) point is M31 globular cluster B29 (see text for details). 



21 



0.4 



o 

-+-> 
cd 

o 

I 1 

o 

cd 
-a 
a 





: Fe52 „ 




• « .r 
- • • •.'<*• 




— ' •<>*•• 






: , I 1 , I I 


: 



0.15 




[Fe/H] 



Fig. 6. — The relationship between published line [Fe/H] values and the line indices in this study. Metallicity 
sources are as follows: • Huchra et al 1991, o Barmby et al 2000 and □ Bonoli et al 1987. The dotted lines 
represent linear fits to the data, weighted by l/of Fo/H1 . 



22 
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Fig. 7. — The relationship between published line [Fe/H] values and the calibrated WYFFOS metallicities. 
The error bar in the upper left represents the median formal error quoted on the WYFFOS values. A linear 
fit, weighted by the inverse-square of the individual uncertainties, is shown by the dashed line for comparison 
to the unity relation (solid line). The slope of the fit is 0.94 ± 0.02 and the RMS of the fit residuals is 0.24 
dex. 
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0.5 1 1.5 



(B-V) 

Fig. 8. — The color distribution of the M31 globular clusters for which we have spectroscopic metallicities 
(shaded) as compared to the overall color distribution of the GCS (solid line). There is no sign of color bias 
in the metallicity sample as compared to the overall M31 cluster population with available (B — V) color 
information (from Barmby et al. 2000). 
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Fig. 9. — The metallicity histogram for the M31 cluster system (top) and the Milky Way GCS (bottom) for 
comparison. The shaded area in the top plot represents the M31 WYFFOS data. 
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-2 -1 

[Fe/H] 

Fig. 10. — Adaptive kernel fits to the M31 metal-poor and metal-rich GC populations (dashed lines), 
separated according to the KMM mixture modeling results. The sum of the individual fits is shown (solid 
line) along with its 90% confidence interval (dotted lines). 
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Fig. 11. — The spatial distributions of the metal-rich and metal-poor cluster populations. At an M31 distance 
of 770 kpc, 5' corresponds to just over 1 kpc. 
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Fig. 12. — The radial distribution of M31 globular cluster populations with metallicities. 
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Fig. 13. — Apparent V magnitude of the clusters as a function of galactocentric radius, showing the metal-rich 
clusters (filled circles) and the metal-poor clusters (open triangles). 
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Fig. 14. — The observed luminosity functions of the metal-rich and metal-poor clusters within the inner 15' 
of the galaxy are represented by the solid lines. The outer cluster LFs, scaled as described in the text, are 
shown by the dotted lines. 
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Fig. 15. — [Fe/H] as a function of projected radius for the full sample and metallicity populations of the 
M31 GCS. A sliding-bin fit was used to determine the mean and RMS at each data point (shown by the 
solid and dotted curves, respectively). 
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Fig. 16. — The velocity histogram for the M31 cluster system. The solid-line histogram is that for the 
available sample of 321 M31 globular cluster velocities, while the shaded area shows our contribution of 191 
WYFFOS velocities to the overall sample. The M31 systemic velocity of —300 ± 4 kms -1 is shown by the 
dashed line. 
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Fig. 17. — Radial velocity of the M31 globular clusters versus projected radius along the minor axis (left) 
and the major axis (right). Culled robust fits to the data are shown by the solid lines. 
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Fig. 18. — The projected rotation curve of the M31 GCS and its 90% confidence limits are shown by the 
solid and dotted lines, respectively (see the text for details). The dashed line represents the stellar rotation 
curve from Rubin & Ford (1970) for comparison. 
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Fig. 19. — The rotation curve of the M31 GCS based on robust biweight estimates of the mean ("location") 
in radial bins of 10 arcminutes for the full dataset and separate metallicity populations. The error bars 
represent the 90% confidence interval. The solid line shows the best fit to the data assuming solid-body 
rotation in the inner 25 arcminutes and flattening beyond. The derived rotation amplitudes are provided in 
the lower-left corner of each plot. 
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Fig. 20. — The projected position angle of the rotation axis (solid line) and its 90% confidence bands (dotted 
lines) as a function of radius for the full sample, as well as for the metal-rich and metal-poor populations. 
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